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DETERMIlTATIOBr 07 STRESSES AIJD DS^ORHATIOlfS 
OE AIRCRAET PROPELLERS* 
By Eriedrich Seewald 



The propeller is probatly one of the moGt highly 
stressed parts of an aircraft. It is therefore surpris- 
ing that the strength of aircraft propellers ha.s hitherto 
received so little attention. While our knowledge of the 
strength of airplane wings has "been considerably increased 
in recent years, relatively little is yet known regard- 
ing the stresses undergone "by aircraft propellers. It ap- 
pears important therefore to give this subject some atten- 
tion. The incentive to the following investigation was 
furnished by a series of accidents caused by propeller 
failures in flight, A large proportion of these failures 
occurred with a type of propeller which had a peculiar 
curving shape. This propeller could not be made to fail 
on the test stand, even with considerable overloading, 
although a blade root of one of the propellers used was 
already damaged. 

At first thought this fact appears rather remarka- 
ble, since the forces generated in stand tests at in- 
creased revolution speeds are much greater than in flight. 
In particular, the thrust forces on the stand r.re a mul- 
tiple of those in the air. If, however, the action of the 
forces on a propeller blade is more closely considered, it 
is immediately recognized that the shape of the blade has 
a preponderant effect on the bending moments. The cen- 
trifugal force acts lengthwise of the blade, while the 
aerodynamic force acts perpendicularly to it. Under pres- 
ent conditions of propeller operation, the centrifugal 
force is one hundred or more times as great as tlie aero- 
dynamic force. If the blade is not quite straight, the 
points of application of the centrifugal forces to the in- 
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dividual "blade elements do not lie in a straight line and 
consequently generate tending moments which tend to 
stretch the "blade. Since the centrifugal force is prepon- 
dera.ntly great in comparison with the other forces, even 
slight deviations from the straight line are important. 
Of course it does not matter whether the deflections were 
ori{<inally present or were produced elastically hy the 
loading. This is a well'-known fact, and patents have al»- 
ready "boon :-;ranted for a propeller with "blades intended 
to 00 so shaped that the centrifugal and aerodynamic mo- 
ments will just offset one another. The "blade axis could 
then "be regarded as a flexible line on which all the 
forces act instead of on a solid "body. G-reat caution 
must "be exercised, however, with such propellers, as il- 
lustrated "by the above-mentioned accident s« . It is not 
sufficient for one to have only a qualitative knowledge 
of the effect of the forces, "but he must also "be alDle to 
calculate the magnitude of the "bending moments with con- 
sideration of the shape, and elastic properties. The on- 
ly scientific pu'blication known to me on the strength in- 
vestigation of aircraft propellers is the work ty Heiss- 
ner, "Ueoer die Biegungsbeanspruchung von Luf t schr auben 
und die entlastende TiTirkung der Zontrif ugalkr af t , Tech- 
nische 3orichte, Vol. II, rlo. 2> 1917, page 315. In this 
work the elastic deflections and the centrifugal moments 
due to these deflections aro determined for a propeller 
which is conceived from the first as a straight untwisted 
"bar, for which a definite di str Ibuti on of the cross sec- 
tions and inertia moments, is assumed. Proceeding from 
similar assumptions, Wilhelm Hoff developed a graphic 
method for determining tho strength of propellers which 
he explained in his lectures at tho Borlin Polytechnic 
Institiito. This work has not yet "been pu"blished (1927). 
Based on these works, C. Jansen undertook an investigation 
on how tho strength calculation of a propeller is affect- 
ed Idy twisting tho blade sections (D,V,L# Report Ho. 50, 
Zoitschrift fur Plugtochnik und Mo t or lui t schif f ahr t , 1925, 
page 07)-, The assumptions concerning the shapo of the 
blade axis and the distribution of tho inertia forces 
along the blade cro essentially the same as in Rcissnor^s 
work. Since, however, for the above-mentioned reasons, 
even slight deviations from tho preliminary assumptions 
entail considerable changes in tho bending moments, it 
is obviously important, in propeller investigations, to 
introduce tho true form into tho calculation, especially 
because, in most propellers and for constructional roa- 
sons, the original blade axis deviates more from a straight 
line than is ascribable to the elastic deflection. 
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In what follows, we must therefore regard the air- 
craft-propeller blade as an originally "bent and twisted 
bar, which is still further elastically deformed under an 
added load. The accurate solution of this problem en- 
counters difficulties out of all proportion to the attain- 
able results. The whole problem, however, can be greatly- 
simplified . by remembering that only the excessive magni- 
tude of the centrifugal force, which acts approximately 
parallel to the longitudinal aris of the blade, malces it 
at all necessary to consider the very slight original cur- 
vature. Hence, in determining the conditions of eqr.ilib- 
rium, we will consider the total curvature (embracing that 
elastically produced and that originally present) only in 
so far as it involves the resolution of the static forces, 
which are so large in comparison with the others that even 
a relatively small component has some effect. If the cen- 
trifugal force, which is approximately parallel to the 
axis of the blade, is resolved into a force lengthwise of 
the blade and another force perpendicular to the blade, 
the latter component is only a .small fraction of the cen- 
trifugal force. Since, however, the other forces acting 
in this direction are very small in comparison with the 
centrifugal force, even this small com\)onent plays apart. 
Conversely, it would be useless to consider any small 
com-ponent of the already small -aerodynami c force, if it 
should fall in the direction of the centrifugal force. 
Moreover, since the deviations of the blade axis from a 
straight line are very small (in extreme cases not over 
1/30 of the blade length), the blade may be treated, with 
respect to the forces perpendicular or nearly perpen- 
dicular to the blade axis, as if the blade axis were 
straight , 

Even in calculating the internal stresses, the cur- 
vature may be disregarded, and the formulas of elementary 
mechanics (for example, a = My/j) may be considered val- 
id. The resulting error is negligible since the radii of 
curvature are multiples of the cr 0 s s- sect i onal dimensions, 
Kote that we are here dealing with cro'-3S sections which 
are asymmetric with respect to the axis passing throiigh 
the center of gravity and that, eveii in this point, the 
assumptions of elementary mechanics are therefore not 
strictly fulfilled. However, since the deviations from 
symmetry are not excessive, we may expect to obtain suf- 
ficiently accurate results. 
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COITDITIOITS OF EQUILIBRIUM 



The tlade axis is the lino passing through the cen- 
ters of gravity of the individual sections. It may be 
any curve in space, but deviates only a lit tie from a 
'Straight line. We shall establish this curve by the fol- 
loAving system of coordinates^ The z axis is placed in 
the radial direction so as to conform as closely as possi- 
ble to the blade axis, the y axis in the direction of 
flight and the z axis perpondictilar to both. The pos- 
•itive directions are so chosen as to form a right-hand 
system. The centrifugal force then acts very closely in 
the direction of the x axis. The centrifugal force de- 
veloped by a blade olomont of length As is designated 
by AC« The aerodynamic force is assumed to act in tho 
yz plane and is resolved into the components Py and p^. 
Wo then have:- 

R g R 

J p^ dx = — and / p« dx x = — 
0 ^ ^ 0 



a 



in which S denotes the thrust, M the torque of the en- 
gine and a the number of blades. 

Two projections of the blade axis are shown in !Fig- 
ures 1 and 2. The manner of resolution then follows di- 
rectly from the above statements, especially the justifi- 
cation for designating the forces in the directions y and 
z as transverse forces. The component in the direction 
of the blade axis is designated by AH and the components 
in the directions of the y and z axes by AQ,y and AQ^g > 
respectively. The magnitudes of these forces are: 

\ c 

AN = — ~ — r- = AC 
cos(sx) 

A9,y = Py Ax - AC GOS(Gy) 

AQ2 = Pz Ax - AC cos(sz). 

The cross section of the blade at the point x = x^ is 
acted on by the normal force 

H 

K = S AC (1) 
and by the transverse forces: 
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R R 
: Qy .- S Py Ax - S AC.cos(sy) (2) 

■ , , . Xi Xi 

R R ^ 

-Qz = S Pz Ax - E AC cos(sz) (3) 
. ■ . -^1 . ^1 

The components of. the bending moment are denoted by My 
and . The bending moment My tends to produce rota-- 
tion about the y axis. It is positive when the result- 
ing curvature of the blade axis is concave in the positive 
z direction. Corresponding statements apply to Mz . The 
components of the bending moment are; 

R • R " . 

% = / py(x - Xi) dX ^ / (y yj dx (4) 

Xi X 1 

R R ■ 

Mz = / Pz(x - dx ^ / . (z - zi) dx (5) 

Xi Xi 

The torsional moment is due partly to the fact tha.t the 
aerodynamic force acting on a blade element does not pass 
through the center of gravity. Shifting the aerodynamic 
force toward the center of gravity would require a corre- 
sponding moment which may be designated by Mod» The 
torsional moment is also due in part to the curvature of 
the blade axis. Hence, 

R dM^ ^ 
Md^ / "TZ'Ccos (sx)-cos (si x)] dx-cos(six) / p (z-Zi)dx + 

xi xi y 

+cos(six) / Pj- (y-yi )dx+cos( Si y) / [p (x-Xi)^ T5(y--y, )] dx - 
Xi ' Xi AX 

R 

-cos(siz) / [Pz(x-x^)^ AC (z~zi)] dx (6) 
xi 

The last two expressions represent the components of the 
bending moments My and Mz , in the direction of the blade 
axis. Since the torque is small in comparison with the 
bending moment, these componentsi must be talien into con- 
sideration, while the torque components in the direction 
of the bending moment may bo disregarded. 

The bending moments are the most important in the 
strength investigation and must be resolved in the direc- 
tion of the principal cro ss- sectional axes. The determi- 
nation of the principal axes is quite troublesome. For 
practical purpo se s i t " suf f i ce s to take one direction 'oar- 
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allel to the "blade chord and the other perpendicular to it 
The comparison of the computed principal directions for 
several hlade sections with the assumed ones yielded a dif 
ference of only about two degrees, which is of no practi- 
cal importance. There is naturally no ohjection to using 
the actual principal axes in the further consideration of 
the' subject. They will be designated by I and II, as in 
Figure 3, The angle formed by the projection of axis I 
with the z axis will be designated by a. This angle 
is called the pitch angle and is defined by the expression 
tan a = H/Sirr, where H is the pitch of the propeller. 
On resolving My ^'^'^^ in the direction of "the princi- 

P8.1 axes I and il , we obtain: 

Mj = My cos a + liz sin a 



Mjj = - My sin a + li^ a. 

The elastic curvatures thus produced are designated by 
kj and kjj corresponding to the inertia moments Jj 
and Jjif and we then have: . 

M 1 

kr = = (i;iy cos a + Mz sin a) (?) 

^ E Jj E Jj 

kiT — ^ii- = — — (- My sin a + Mz cos a) (8) 



If S- denotes the angle of elastic distortion of the blade 

and the torsional strength, the torque will be given 

by the expression 

= (9) 

dx Gr Jd 



The c-arva,tures k and the torque here represent sim- 
ple elastic deformations, which must be expressed by the 
coordinates of the blade. Kence the coordina-tes must be 
divided, on the one ha.nd, into those representing the o~ 
riginal shape and, on the other hand, into those repre- 
senting the elastic deformations. The shape of the un- 
stressed blade axis is denoted bjr the coordinates 

^o» J'o^ ^0 ^0 
and the corresponding elastic deformations by I, T], 
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^ and 1^. The variation in the length of the blade is es- 
sentially represented by i and is so small that it may 
"be disregarded. 

The final coordinates of a point on the blade axis in 
the stressed condition are therefore 

x = Xo,y = yo+^tZ = 2o + ^ and a = ao + n^. 

The elastic curvatures kj and kji, expressed by these 
coordinates, are: 

k = _^ COS ao + ^) + ~ sin(ao + = —7- 
dx^ cLx^ E Jj 



^11= - y Bin(ao_+ ^) + ^-^ cos(ao + ^) = ^ • 

In equations (7) and (8), the quantities Mj and Un are 
expressed by My and Mz . By introducing these va.lues, 
we obtain 

ki = iln cos(ao + ^) + 7^ sin(ao + ^) = 
dx2 dx^ 

My cos(ao + S-) . I'^z sin(ao + iJ) 



+ .-I^ (10) 

E Ji E Ji 

- ^ sin(ao + 1?) + 1^ cos(ao + ^) = ~ 
dx^ dx'^ 

My. sin(ao + _^ cos(ao + ^) ^^^^ 



To these we must also add the equation: 

dx G Jd 



(12) 



For convenience the value of is not written out in 

'equation (S) . The moments My and Mg are expressed in 
equations (4) and (5) by the coordinates of the blade a.xis 
and by the forces acting on the blade. If, on the one 
hand, the coordinates y and z in the equations are di- 
vided, as mentioned above, into the components y^ and Zq 
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which determine the original form and', on the other hand, 
into the components r, and ^ which denote the elastic 
deformations, equations (4) and (5) then become 



= / Py(x - Xi) dX - 



E 

l-^Z = / Pz - ^l) 



xi 

R 



- / 4^ (ri - TiO dx 



R 



Xl 



Ax 



(13) 



■ / f - to dx 



(14) 



These values miist he introduced into equations (lO) to 
(12). in order to obtain the conditional equations for the 
elastic deformations in ter.ns of the blade form and ox the 
acting forces. In equations (13) and (14) the first t^To 
integrals represent the moments which would be produced b^ 
the forces acting on the unbent blade. These can be cal-" 
culated directly. The last integrals in the equations for 
My and liz contain the change in the moments due to the 
change in the shape of the whole blade under stress and 
the consequent variation in the points of application of 
the forces. The nature and magnitude of these deforma- 
tions can be represented by exponeiitial series a,nd write 
the e:Kpr e s si ons : 

T) = ao + a^ X + a^x^ + a3X^ + 

^ = bo + bix + bsx^ + bax^ + 

^ - Cq + CiX + Csic^ + CaX^ + 



Wo calculate with these expressions as if they were l^nown 
quantities, lea,ving the constants a, b, and c to be sub- 



sequently determined so that 
be represented as accurately 
If we introduce these series 
we can calculate all unknown 
b, and c. If we compute numerically 
these equations for any values of x 



the actual deformations will 
as possible by these series, 
into equations (lO) to (14) , 
quantities as functions of a, 
the quantities in 

we obtain, for each 



value 



of 
and 



three equations in which only the quantities 
occur as unknowns. If, for example, we wish 
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to consider the first three terms of each series, we must 
compiite equations (10) to (14) for three values of x and 
ottain -nine -equations for the nine uniinowns. In this 
form the solution would still "be very troulDlesome , 
the unlinowns occur in every equation. This 
remedied as follows. Tie first assume that 



thu 

general 
because all 
difficulty is 



in the angular 



1^ = 0, leaving only CL^^ 
the original twist. Furthermore, we put 
justified "by the fact that, in propeller 
ertia moment Jjj is always much larger 
will appear later that "both S- ^i'^^ kjj 
small that they can he put equal to zero 
proximation. Equation (ll) then "becomes 



f unct i ons , i.e., 

kji = 0. This is 
sections, the in- 
than Jj. It 
are actually so 
in the first ap- 



-II 



sin an + 



dx'^ 



cos 



(15) 



On multiplying this equation "by sin ao ^nd 
it from equation (lO) multiplied by cos ao » 
likewise i5 = 0, we obtain... 



subtracting 
in which 
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(My cos^ ao + ^^z sin a cos ao) 



E Jj 



(16) 



The only other quantity now appearing on the left side is 
the second differential coefficient of T]. The right side, 
however, contains, in addition to T), as follows from 
equation (5), also ^ in- . By assuming kjj = 0, w.e 
can also- express ^ as a function of T] and. therefore by 
a. Equation (15) then yields 



= dill tan 
dx^ dx"' 



a. 



As. already explained, a is the pitch angle and conse- 
quently, 

H 



t an a = 



2 IT X 



For most propellers the pitch H is approximately constant 
over the whole blade. Where considera.ble deviations occur, 
tan a should be expressed approximately by a correspond- 
ing function of " x and the procedure should then be exact- 
ly as follows: 



dx2 



dx^ 2 rr X 



(17) 



10 



IT.A.C.A. Technical Memorandum lid • '670 



For r\ we developed the expression 

T] .= + a^x + a^x^ + ^3^^ + 

We can immediately put ao = ai = 0, "because the "blade 
may te regarded as fixed at the hub), so that T] and dr|/ 
dx 'must vanish for x = 0. Hence, 

=^ 2 a2 + 6 agX + 12 a^x^ +.20 asX^ + 

dx^ • 

By introducing this value into equation (17), we obtain 
= (2 ap + 5 a^x + 12 a.x^ + 20 agx^ -f . . , ) — . 

dX^ ^ ^ 2 TT X 

^ can now "be computed by integration. 

= (2 as In + S a3X 4- 6 a^x^ + 

dx 

+ f ^ ^ ^ + 

d^/dx must equal zero for x = 0, which is possible only 
when a^ = Ci = 0. 3y another integration, we obtain 

^ = (3 aaX^ + 2 a^x^ + I asX^ + . . . . . •) ~ (I7a) 

How all the quantities in equation (16) can be expressed by 
r\. From this equation we then calculate an approximate 
value for r\ and simultaneously also for ^. We can then 
compute all the bending moments. At first these values 
are only approximate. In the numerical example, however, 
it will be found that the assumptions kn 0 and = 0 
are so accurately fulfilled, that the values calculated on 
the basis of this assumption may be regarded as final. 
Should this be otherwise in particular instances, the val- 
ues may be corrected by introducing the calculated approx- 
imate values S- and kjj and making the calculation again. 

Hence the problem is to determine the constants a 
of the series 

= % x3 ■ + ai x^ + a5 x^ + 

in such a way that equation (16) is satisfied for as many 
values of x as there are terms in the series for r\. 
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We "break off the series for r\ with, the third term and 
find, as the first condition, that all moments "become zero 
at the blade tip, including also the curvature d^ri/dx^. 
The other two conditions are determined as follows. If 
a,n assumed function is correct, and the "bending moments on 
the right side of equation (16) are calculated, these "bend 
ing moments necessitate certain definite deflections, ■ If 
these are calculated, they must equal the assumed vali-'es 
of at every point. This condition cannot he reajized 

at all points, hut at as many points as there are toicps 
in the series. Since we have not yet had over two con-- 
stants, wo can satisfy this condition for only two points. 
vlo arbitrarily stipulate that the deflection at the blade 
tip and for x = 0.6 R shall be the correct one. We have 
then represented the deflection by an approximate function 
which agrees with the actual bonding lino under the iollow 
ing conditions. At the hub it is fixed, while it has the 
correct deflection for x = 0.5, At the blade tip it has 
the correct deflection and the correct curvature. It may 
be assumed that such a function, which satisfies five 
boundary conditions, will not deviate greatly from the cor 
rect curve. As conditions for the last two constants, it 
was stipulated that the deflections and not the curvature 
should agree at two points. If the curvature had been 
chosen, then two conditions for the constants would have 
been obtained. The result would still have been very in-* 
accurate, principally due to the fact that, on the right 
side, only the inertia moments were introduced at the two 
places. As a result of the double integration of m/eJ, 
however, the magnitudes of the inertia moments and cross 
sections were indirectly taken into consideration at all 
points by means of the deflections. 



NUMERICAL EXAMPLE 



The process of calculation will be illustrated by an 
example, taking the propeller described in the introdtic- 
tion. It will be found that the stresses in flight, espe- 
cially at high speeds and therefore at low thrust, are 
considerably greater than at the same revolution speed on 
the stand. The propeller used in the tests was of this 
type with a diameter D = 3.45 m (8.04 ft.) a,nd a pitch 
H = 1.15 m (3.77 ft.). It was tested on the stand at 1450 
r.p.m» and yielded a thrust of 325 kg (716.5 lb.) with an 
engine torque of 47 ^m-kg (340 lb. -ft.),. The stresses were 
determined under these conditions. The propeller was run 
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at a speed of 1700 r.p.m. on the stand without injury. 

The pr opeller-"blade sections were 10 cm (3.94 in.) 
apart. In Figure 4 the cross-sectional area is plotted 
against the radius. The e.g. of each measured section was 
found and the coordinates and of the "blade axis 

were determined. The "blade-axis curves are shown in Jig-- 
ure 5. Figures 6a and 6d show, respectively, the sine and 
cosine of the angle made "by the "blade-section chord with 
the plane of the propeller. Figures 7 and 8 show, re- 
spectively, the inertia moments Jj and Jjj at each point 
of the "bla^de (Jji on ten times the scale of - Jj) . All 
these quantities are plotted against the nondimensional 
abscissa I = x/r, where x represents the distance of 
a cross section from the center of the hub. 

Equation (16), when the above-given series and the 
substitution i = x/e are introduced for r\ and I with 
consideration of equations -(13) and (14), takes the form 




. / AC (y -y 00.^1 [H^^ a -^1^) 



+R^a^(r-e,^)+H'a3(e'-e,^) ] 



+ si ng^..^ f / ^ Pz(e-ex)cie-/^ (zo-zoi)de} 
- sing^cosg / M [0.55 a^(f + 

E Ji A^ . 

+0.366 R^a4(e^-e,'')+0.305 RasH^-l^^) 1^^ (la) 
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In order to determine the deflection itself, it is neces- 
sary to integrate tv/ice, wheroty it must "be remembered 
that a and J are functions of i. All the terms on 
the right side can n.ow "be .comput.ed, Hpwever., the integra- 
tions cannot all be computed since a, J and the centrif- 
ugal for.ce at Qp.cn po.int ar.e not .all giyen in. terms, of i , 
but take a, course which can be determined only by measure- 
ments. Ail these expressions are therefore obtained by 
summation. The distribution of the aerodynamic forces is 
approximately known, For the strength cal'cul at i on , it 
suffices to approximate them by a simple function. For 
example, the component forces which produce the thrust, 
may be assumed to have a distribution which is approxi- 
mately represented by the function 



p^ = k i\ri - 

The constant k must be so determined tha,t the sum over 
all the blades will give their total thrust. It m^i.st 
therefore be 

Q \ d. 

whero a is the fluid coofficiont, whonco 



k = 



16 S 
a TT R 



The 'bonding moment duo to the thrust is theroforo easily 
calculated. At the point i - ix i* is 



1=1 



C = ^i ^ TT n ^=--5^ 



After integrating and again replacing t>y ^> ob- 

t ain 



a TT 



^~..-L_ (6 + 7 + IS) - 
120 



- it (I- arc Sine); 
Hence the "moment exerted by the thrtist is 



a TT. 



14 



N.A.C.A. Teclmical Llemor andiim Ho. 670 



in which the oracket^d portion of the preceding formula is 
replaced "by f(|)» .This function has the values: 

i = 0 0.1 0.2 0.3 0.6 0.5 0.6 0.7 

f (0=0.133 0^.114 0.094 0.075 0.057 0.040 0.0255 0.0139 

i = 0.8 0.9 1.0 

f(i) = 0.0055 0.0011 0 

The aerodynamic component acting in the plane of the pro-- 
poller disk is represented accurately enough hy the func-* 
tion 

We again determine kj , so that the total moment exerted 
"by P2 is equal to the measured or calculated engine 
torque. 



^ = f ^ i H = K ^ f e' y 1 - i" de - k, 

a o 0 . 16 



Hence 

k 



Mm 16 



^ TT a 



Pz = 



Mm 16 ^ ^ j,2 



IT a 



Then the "bending moment due to this aerodynamic component 
at any point is 



Mm 15 
TT a 



(2- e - 1) 



- — arc sin i + 
8 



TT 



For ahhreviation the oracketed portion is replaced "by 
q)(^). It has the following values: 
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i = 0 0.1 0.2 0.3 0.4 ■ 0.5 0.6 

cp(O=0.1962 0.163 0.1312 0.106 0.073 0.0502 0.031 
i = 0.7 0.8 0.9 1.0 . 

cp (O = •0*C)167 0.0044 0.0014 0 

7ith these the "bending moments due to the aerodynamic, 
forces are determined. Their curves are plotted in Fig- 
ures 9 and 10. We now have to determine the moments due 
to *the centrifugal force and first tliose' due' to the orig- 
inal curvature. These "bending moments are represented 
"by the expressions: 

4^ (yo - yoi)Ae and 4f (^o' - ^oi)H • 

The quantities yo. - 'Vi ^^^^ .^o ~ ^'"^^ '^^ read from 

figure 5 for any value of 

. If we imagine the "blade divided perpendicularly to the 
i a.xis into' disks of the thickness 1 cm = , then the 
centrifugal force per unit length is AC/RAE- 

= ^ 1 R i (j^^ 

R A r g 



The. specific gravity of the wood is assumed to "be 0.8 or 
800 kg/m^ (28 or 28250 cu.ft.). At" 1450 r.p.m. Oj^ = 
33 000. In order to oh tain we must therefore mul- 

tiply -the curve, which indicates the cro ss- secti onal area, 
"by the factor 



S 



I = 8 00 X 2300 0 X 1.9 25 ^ ^ 229 | 



9.81 X 100 



unit length of the blade is plot- 
calculation, the "blade is di- 
planes perpendicular to the i 
0..1. I.n order to determine the 
such, portion , we must determine 



This centrifugal force per 
ted in Figure 11. For the 
vided into sections by ten 
axis at interva,ls of A^ = 
centrifugal force of every ., ^. . 

the* surface area of. the. curve. AC/RA^ in Figure 11 between 
the limits ^ "and.. I + At. We may conceive of such a por- 
tion with. s'uffici"ent accuracy ,as a .trapezoid, 



Figtire 11. The point of a.pplic at i on of this 



as snown m 
centrifugal 
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force to the hlade then lies at the value ©f I which cor- 
responds to the e.g.* of this surface strip. In what fol- 
lows, however, it is more cgnvenient to resolve this force 
in such a way 'that all the components are applied at 
i = 0,1, 0.2 , etc. This resolution can oe made "by apply- 
ing half of the area of the hatched rectangle in figure 
11 to each of the. two points I. and i ^ 1$ and 2/3 of 
the area of the triangle to one point and the remaining 
1/3 to the other point. The centrifugal forces for all 
the "blade elements were thus determined. 

i = 0 0.1 0.2 0.3 0*4 0.5 0.5 0.7 0.8 0.9 1.0 

56 358 769 1081 1207 1225 1216 1140 949 602 110 

For determining, at the point ^, the "bending moment due 
to these forces, we have the expression: 



if (jTo - y..) Ae or Jz'^ If - - 



^n 



01 



) Ae 



The simplest way to solve such an expression id "by using 
the formula, 

Mn = Mn-i + A yon ^ A e ^ • 



Ayon Azon are to "be taken directly from Figure 5, 

where the original "blade form is plotted. The moments 



Mq^^ thus obtained are plotted in Figures 12 



and 13. We are now in position to solve a portion of the 
right side of equation (18), namely, the expression: 



The "bracketed expressions have alread^r been solved and 
represent the bending moments- due to the action of the • 
aerociynamic and centrif^igal forces on the unbent blade. 
The angle a and the inertia moment J are measured at 
each point. E = 100000 kg/cm^ (1420000 Ib./sq.in,) was 
adopted as the elasticity modulus of wood. This whole ex- 
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pression has to he integrated tv^ico, in order to ohtain 
the deflection, which would he produced hy the "bracketed 
moments. A double integration is equiva^lent, however, to 
the production of a moment according to Mohr * s theorem* 
The quantity represented by the above expression may there- 
fore be taken as the load applied to a blade, from which 
the bending moment can be calculated. The evaluation can 
be made exactly as in the above , determinati on of the bend- 
ing moments. The result is plotted in ^'igure 14. The 
quantity, calculated as a, bending moment, is designated 
by and represents the deflection produced by the 

above--calculated moments, wlxen the moments remain constant 
during the deformation. The two still-unsolved expressions 
in equation (18) 

i^) +0.366 a4(f - ei"") + 

+0.305 R^ a5(f - ^i") ] d ^ 



and / |§ [0.55 R^ ^-sif - 



represent only the moments which are first produced by the 
elastic deformation and -which overlap the above-calcula.t ed 
momen t s , 

We will now continue with the terms containing the 
constant as. On extracting these terms we obtain* the two 
expressi ons : 

C.55 a3 r2 / 11 _ e,^) de. 

Only known quantities stand under the integration sign. 
These can 'therefore be calculated by summation for any de-- 
sired value of the calculation being made the same 

as for the previous' quantities,' as 

II (yo yoi) ' 

Both integrals represent functions of ^ which can be 
directly computed and which we will represent by cp(^) 
and (i) . a3 'pH) and 0.55 as (i) then represent the 
bending moments produced when the blade is bent into the 
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form T| = a3 ^ . By introducing- tlie se quantities into 
equation (18), we obtain 



9'3 



^osa ,r:(l) + 0.55 illL^^93^ ^(l) 
E J J E J J 



as the component of the oend produced by these moments, 
Again considering the bracketed expresa on as the load and 
having determined the bending moment, the latter is equlv-- 
alent in magnitudie to the deflection produced by the mo- 
ments 9(0 simd 0.55 i/ii) in the y direction. This 
deflection is designated by 1^3. The procedure is the 
sr.mo with the terms. . and as, and the corresponding 

deflections are represented by r|4 and ris. 

All these deflections are plotted in Figure 15 on' a 
magnified scale of 1000:1, The total deflection is then 
equal to the sum of all the individual deflections, so 
that the equation must be satisfied for every^ point ^x. 

T] = a3 £3 ^^3 4, 4- a^ - 

As r.lready explained, we put l-^ on the left side, once 
ii =, D.S and once ,ii = 1> while we introduce the calcu- 
l.a.ted deflections . T\q , 7^3, r]^ and r\s on the right side. 
Thus we obtain two equations for the constant" a. For 
the third equation, as already mentioned, we choose the 
condition that, at the blade tip (that is, for - 1) $ 

all moments equal zero, so that d^ri/d^^ miist vanish. 
Thus we obtain the three equations: 

1.225X1 + 12 a4 1.225^ + 20 1.225^ = 0 

100G(a3 1.225^X0.6^+a4 1 .225^ XO . 6^ H-a^ 1.225^X0.6^) = 

=l,473-a3 250.16^a4 254,7-a5 247.6 

1000(a3 1 .225^+a4^ 1.225^4-ar3 1 .225^) - 13.126 - 

-a3 854,9-.a4 924-a5 1062.6. 

These equations yield, for the a constants the values . 
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as = -0.017555 



a^ = +0 .036033 



as = -0.01414 



Hence tlie total c'eflection is 



n = a.z'R^i^ + a.-.R^I^ + a5R°|^ = 
+ 0.03603 X 1.225"* 
Accordine to equation (l7a) 



0.01755 X 1 .225^ P + 
i'^ - 0.01414 X 1.225^ 



t = (3 + 2 a^ |^ 



The quantities r\ and ^ are plotted in figures 16 and 17, 
from which it is obvious that the calculated deflection I 
is very small as compared with the deflection in the y 
direction. The shape of the propeller, after its defor- 
mation, is now known (fig. 18), and all moments can be 
calculated by equations (4), (5), and (5). The bending 
moments for as = 1, a^ = 1, and as = 1 have already 
been calculated, so that we only need to multiply them by 
the computed values of a. The moments Mj and Mjj can 

then be readily calculated from the individual expres- 
sions, Likev/ise the stresses can be directly calculated. 
The results are plotted in figures 19 to 22. The bending 
stresses due to Mg v;ere also calculated, though they 
are so small, due to the magnitude of Jg, that they 
play no part, all the more since the maximum stresses due 
to Mjj occur in the fibers v/here the stresses due to 

Mi are almost zero. Since all the quantities are now 
known, the calculation of the shearing stresses can like- 
wise be made as for ordinary beams. These stresses also 
are not very important. The torsional moment can be read- 
ily calculated from equation (6). It is so small, however, 
that it is of hardly any importance. It is plotted in 
figure 23. All these quantities are now calculated on the 
assumption that the curvature in the direction of the 
greater inertia moment and the elastic distortion of the 
blade are both negligibly small. In order to determine 
how far these conditions are satisfied, v/e v/ill now cal- 
culate them. Since the magnitude of a distortion for any 
given cross section is very difficult to calculate, every 
cross section is replaced by an inscribed rectangle. The 
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torsional rigidity of such a "blade can then he determined 
according to St. Venant and is surely less than that of 
the actual hlade. Likewise the modulus of shear G- is 
assumed very low, G = 8000 kg/cm^ (113800 Ih./sq.in.). 
Under these assumptions, we have, at the propeller tip, a 
twist of 75 = 0.00661 or abo.ut 0.3°. This angle, howev- 
er, is within the degree of experimental accuracy. The 
torsional an^'le night he somewhat greater for other "blade 
forms. Nevertheless, it would hardly reach values which 
would affect the strength relations. 

In order to test the other assumption kji = 0, we 
will compute kji and compare it with kj. These quan- 
tities are calculated for a few values of \ "by forming 
the expressions Ki/EJj and Mji/EJjj. 



i 






= 0.2 


-1667.0 


+ 29.2 


= 0.4 


+ 2620 


-17.4 


= 0.6 


+ 5140 


-36. 8 


= 0,8 


+ 4430 


-20.0 



It is obvious that kjj is much smaller than kj. At the 

r.ost unfavorable point (| = 0.4), kjj = 0.0175 kj. It 

seems perfect ly • permiss ihle , however, to disregard such 
a small quantity. This proves that the assumptions kg = 
0 and = 0 do not greatly impair the value of the fi- 
nal result, so that we may regard the calculated static 
data as correct . . * 

It is obvious from the course of the bending moments 
and stresses that the m.oments due to the centrifugal 
forces preponderate, even at the great thrust on the stand. 
The blade is bent backward at the root even at groat 
thrust. It is therefore of interest to know the bending 
moments when the thrust is small. It is therefore assumed 
that the propeller runs at the same speed as above (U = 
1450 r.p.m.) when it exerts no thrust at all. This is 
a condition which can easily occur in flying at a low an- 
gle of attack and with the engine throttled. If all the 
aerodynamic forces are put equal to zero, the calculation 
is greatly simplified, since all the data required to 
solve the equations have been previously determined. It 
is only necessary to put all expressions involving aero- 
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dynamic forces ^qual to zero, Tlie calculation, is tiien 
made exactly as "before. Since it affords nothing new of 
interest, only the results will te given. The stresses 
are plotted in Figures 24 and 25 and show that, near the 
blade root, they are considerably greater than in the pre- 
vious condition, where the aerodynamic forces act. This; 
.is also quite natural, for even in the previous loading 
case, where the. aerodynamic forces tended- to bend the 
propeller forward, the opposed "relieving" moments, which 
are produced hy the centrifugal force due to the curved 
shape of the blades, so fa,r exceed the aerodynamic moment, 
that the resulting moment is even greater than the moment 
vmich the aerodynamic force alone would exert on a 
straight blade. In normal flight, however, the aerody- 
namic forces are considerably smaller and they may almost 
vanish when the airplane is .in a glide. The aerodynamic 
force then produces a backward-^bending moment, which is 
much greater than the aerodynamic forces could produce, 
even in the most unfavorable case. The blade tends of 
itself to assume such a position that the effect is min- 
imized. The blade can undergo important deformation, 
however, only in its outer part, where the cross sections 
are small • Hear the hub no considerable deformation can 
occur, unless the stresses are very high.. The calculated 
stresses never suffice, however, to cause certain fail- 
ure, even when they are :dangorously. high. It must be re- 
membered, however, that the conditions may become consid- 
erably more unfavorable ^through some small material de- 
fect, such as may Joe produced by the weathering of a glued 
joint. Further considerable stresses may be produced 
by even small vibrations. The previously mentioned fail- 
ures occurred mostly in ^largo propellers, which were sim- 
ilar in shape, however, to those investigated here^ In 
so far as known, all the fail^ires occurred in flight, i.e., 
at low thrust, as might be concluded from the present in- 
vestigation. The cause of the failures, therefore, was 
probably excessive curvature of the blades. 

Hence great caution must be exercised in trying to 
give a propeller such a shape that the aerodynamic forces 
will be counterbalanced by the centrifugal forces. It 
seems particularly hazardous to let the blade axis near 
the hub project too far from the plane of the propeller 
disk. The effect of a faulty shape is not so detrimental 
in the outer portion of the blade, due to its flexibili- 
ty, which partially remedies the defect. Hear the hub, 
however, due to the great rigidity, no considerable • auto- 
matic adjustm^ent is possible, and the moments are great- 
ly affected by the original design. It is impossible to 
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establish any. general rule for the design. It may, how- 
ever, "be said. that, in the case. of *a "blade having a . 
straight axis, tlie elastic deformation is always such 
that the load is. actually reduced "by the centrifugal force, 
regardless of all other forces acting on the propeller. 
This fact was also demonstrated, "by Reissner. For a curved 
"blade this is not the case. In any event it should be de- 
termined whether the load reduction under . one operating 
condition does not involve a load . increase under another 
operating condition. 

The propeller blade has thus far been considered as 
a thin -plate. This was accurate enough for the quantities 
hitherto calculated. This assumption causes the complete 
disappearance of one factor, which should be considered 
in many cases. This is the torsional moment exerted on 
the twisted blade . by . the longi tudinal force. If, for ex- 
ample, we imagine a metal strip twisted about its longitu- 
dinal axis so as to form a helical surface, and, if this 
strip is then subjected to longitudinal tension, it will 
tend to return to its original flat shape. In order to 
determine these moments approximately, each fiber may be 
considered as a twisted blade,. When this is pulled, it 
will- tend to untwist. The . loiigi tudinal stress in each fi- 
ber is known accurately enough from the previous calcula- 
tion. Due to the twist, the longitudinal direction . con- 
tinually varies along the fiber. There must therefore be 
everywhere a component perpendicular to the fiber, and 
this produces a "backward turning moment. In thick propel- 
lers, such as those of wood, this, moment is unimportant, 
but in very thin blades it may. assume great importance 
and cause considerable distortion, .For the approximate 
calculation of this torsional moment in a thin propeller 
'blade, it is siifficient to take account of the centrifu- 
gal force alone (assumed to be uniform throughout the 
cross section), since, due to its slight rigidity, the 
blade automatically takes such an attitude that the bend- 
ing stresses are small. If the angle through which a sec- 
tion is twisted with respect to the preceding section is 
measured at def ini t e int erval s of 10 cm (3,94 in,), for 
example,' this angle is then a criterion for the bending 
of .the individua.1 fibers. The components of the internal 
forces ,. producing the torsional moment, can then be read- 
ily determined from the change in the direction of the 
fibers. 

Various tests were made in order to determine the 
deformation of the propellers experimentally. Though all 
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these tests vrere of a rather primitive nature, it wa.s 
still possible to verify the c£il dilation with one of the 
arrangements tested. (Figs. 26 and 27.) The whole con- 
tour of a "blade of the above mathematically investigated 
propeller was lined with sparL; go.ps by attaching strips 
of tin foil 2.5 cm (about an inch) long. One end of this 
conductor was connected with a condenser by means of a 
slip ring with brushes, while the other end consisted of 
a point which passed a st at ionary . point once during every 
revolution at a distance of about 2 mm (0,08- in*) i tlie 
latter point being connected with the opposite pole of 
the condenser. The energy was so adjusted that a spark 
could pass only when the revolving point was exactly oppo 
site the stationary point. The condenser was charged by 
an induction coil operated by an alternating current. 

•First,. in order to determine the. shape of the propel 
ler in the unstressed condition, the blade was turned to 
the po si t ion , where the spark could pass. A photographic 
camera was then mounted in the plane of the propeller 
disk in such a way that the axis of the lens was approxi- 
mately perpendicular to the axis pf the propeller blade, 
and a spark was discharged through the spark gaps and 
was photographed. The camera was then left in the sa.me 
position, while the propeller was made to revolve. At 
1450 r.p.m. another spark was photographed on the same 
plate. In order to make sure the camera had not moved, 
spark gaps were also arranged at two stationary points. 
In this way Figure 28 was obtained. The deflection can 
thus be determined with great accuracy. A twisting of 
the blade, which would be indicated by changes in the di s 
tances between the sparks on the loading edge and those 
on the trailing edge while the propeller is revolving, as 
compared with the distances between the same sparks when 
the propeller is at rest, can indeed be detected, but it 
is so small that it cannot be accurately measured. A 
twist of one degree would cause a relative displacement 
of the points of 1.5 to 2.5 mm (0.06 to 0.1 in.) accord- 
ing to the width of the blade at the point invo st iga,t ed. 
The deflection is readily recognized, however, and agrees 
well with the above calculation. At the blade tip the 
measured deflection is 10.6 mm (0.42 in.), and the calcu- 
lated deflection is 9.7 mm (0.38 in.)» The general char-- 
actor of the measured deflection agrees well with the cal 
culated deflection. Even the negative deflection near 
the hub is recognizable. In order to test the accuracy 
of this method of measuring, the sparks were photographed 
with the propeller at rest. The blade was then distorted 
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and the deflection was accaratoly measured. at several 
points. Then another photograph was taken and evaluated. 
It Y7as thus found that deformations of 0.5 mm (0.02 in.) 
could he determined, •' 



SUMHAHY 



A method is descrihed for testing the strength of 
propellers of any shape. It is shown that the shape of 
the propeller greatly affects the stresses, and that great 
caution must he exercised if the relievinr^; effect of the 
•centrifugal force, which exists in every propeller with 
a straight "blade axis, is to "be . increased "by curving the 
"blade axis. 



In a numerical example, the calculation is made for 
a certain revolution speed at which the actual distortion 
is then measured. The results of the calculation agree 
satisfactorily with the experimental results. 



Translation hy Dwight M, I.Iiner, 
National Advisory Committee 
for Aeronautics. 
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Figs.l,2,3,4,5,6a,6'b 
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on xy plane. 




Fig. 3 Principal axes and pitch 
angle of 'blade section. 
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against the radius. 
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Fig. 5 Blade- axis curves. 




Fig. 6a Sine of angle made "by chord 
of "blade section with plane 
of TDroDoller disk. 
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Fig. 6b Cosine of angle made 

by chord of blade sec- 
tion with plane of propeller 
disk. 
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Figs. 7, 6, 9, 10, 11, 13, 13, 14 



800, 

Jl 600 
cm 400 
200 



























t 







0 0.2 0.4 0.6 0.8 1.0 

i 

Fig. 7 Inertia moments, Jj 
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Fig. 8 Inertia moments, Jjj 
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Fig. 10 Z component of tending 

moments due to aero- 
dynamic forces. 
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Fig. 11 Centrifugal force per -onit 
len.^th of iDlade. 
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?igs. 15, 15, 17. 18, IS, 20 



•n 



1000 
800 
500 
400 
200 



- - 








1 








tr- 












1 

1 






1. 

A 


i 

1 






-i 




1 - 

I 

1 










— h" 

! 


'J 

1 — .aK" 






1 1 1 

i 



0 0.4|0.8 1.2 
Fig. 15 



l.C — 

0.8 

0.6 

TijCin 0.4 

0.2 

0 



0 0.2 0.4 0.5 0.8 1.0 



]?ig.l6 Calculated deflection 



in direction of flight. 



.08 
.04 

<p , cm 0 

.04 
.08 



I ■■ ■ 








































































/ 








y 


^ 


, — h--^ 







0 0.2 0.4^ 0.5 0.8 1.0 

?ig.l? Calculated deflection in 
plane Of propeller dislc. 



cm 





^ 






/ 






i 1 


/ 




! i i 



0 



0.2 0.4 0.6 0.8 1.0 

b 



Pig. 13 Blade axis with pro- 
peller at rest and at 
full revolution sneed on stand. 
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Fig. 19 Bending moment, Mr. 



Fig. 20 Stress due to cen- 
trifugal force. 
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?igs. 21 , 22, 23, 24, 25, 28 
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Fig. 21 Total stress on pressure side. ?ig.22 Total stress on ten- 
si oii side. 
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Fig. 23 Torsional moment 
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Fig. 24 Stress on pressure side. 
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Fig. 28 Outline of "blade at rest 

rest and at full speed 
on stand (copied from photo- 
graph for greater clearness). 



Fig. 27 Close-up of Fig. 26 
showing details. 




Fig. 26 Experimental arrangement. 



